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Abstract 

Recent research of blind mole rats of the species complex Nannospalax (superspecies leucodon) identified a small 
and fragmented population of these rodents on both sides of the Hungarian-Serbian border. Cytogenetic 
investigations proved that this population karyologically identical with the Vojvodina blind mole rat described 
earlier as Nannospalax {leucodon) montanosyrmiensis. Based on cytochrome b gene sequences obtained from three 
specimens originating from separate locations, these blind mole rats form a discrete phylogenetic clade which, with 
a difference of about 10%, is well separated from other blind mole rat taxa inhabiting the Carpathian Basin. The 
taxon has only two extant populations that are 150 km apart from each other. The combined occupied area is 
estimated to be less than 10 km 2 , and the total estimated number of individuals is less than 300. These two 
remaining populations are heavily fragmented and under imminent threat by the establishment of tree plantations, 
small-scale and agro-industrial farms and land development. The situation is further aggravated by the fact that 
80% of the individuals inhabit unprotected areas. A study of the landscape history of the wider area surrounding 
one of the populations - based on military maps spanning over the last 200 years - has shown a drastic decrease in 
the extent and quality of potential habitats. Based on our present knowledge, the Vojvodina blind mole rat is one 
of the most seriously threatened, rarest mammal in Europe, the remaining population of which can be wiped out 
within years unless immediate conservation action is taken. 

Keywords: Carpathian basin, Conservation biology, Cytogenetics, Extinction, IUCN categories, Nannospalax 
{leucodon) montanosyrmiensis, Spalacinae 



Background 

Rodents (Order Rodentia) are usually not in the focus of 
conservation biology (Lidicker 1989). Proper evaluation of 
the conservation status of rodent species is further compli- 
cated by the confusion that surrounds almost all levels of 
rodent systematics (Corti 2001). Well-established tax- 
onomy is the base for efficient conservation biology, but 
unclear taxonomical questions can easily result in the neg- 
ligence of certain groups, whose extinction would mean 
the loss of entire evolutionary lineages, drastically decreas- 
ing the overall biodiversity on Earth. The vulnerability of 
this order is demonstrated by the fact that rodent species 
contributed 51-52% to mammalian extinctions in the 
last 500 years (Ceballos and Brown 1995, MacPhee and 
Flemming 1999). Counterintuitively though, conservation 
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initiatives will continue to be biased towards the most 
studied and attractive mammal groups and species (Amori 
and Gippoliti 2000). 

The situation is clearly mirrored in the case of the 
Eurasian blind mole rats (Spalacidae: Spalacinae). These 
small mammals represent a distinct group among rodents 
which is extremely adapted to subterranean life. They have 
cylindrically shaped body with no external ear and a ves- 
tigial tail, and are completely blind spending their entire 
life in their tunnel system built underground (Topachevskii 
1969). Compared to other rodents, the conditions resulting 
from their lifestyle created a decreased morphological vari- 
ability and the species are very similar both externally and 
osteologically (Nevo 2000). Putting aside the lineage of 
large mole rats (genus Spalax) (for taxonomic context and 
nomenclatural details see Topachevskii 1969, Nemeth et al. 
2009, Arslan, Akan and Zima 2011, Hadid et al. 2012), taxa 
belonging to Nannospalax present a long-standing source 
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of dispute and disagreement on their systematics (Savic 
and Nevo 1990, Musser and Carleton 2005). Within the 
latter genus one of the recognised species groups (regarded 
as superspecies) which include a large number of karyolo- 
gically different taxa (for the list of these named forms see 
Savic and Soldatovic 1984) is Nannospalax (superspecies 
leucodon) (Musser and Carleton 2005). Molecular genetic 
investigations of this superspecies are quite limited so far 
both in terms of geographic and taxonomic coverage 
(Krystufek et al. 2011, Kandemir et al. 2012, Hadid et al. 
2012) and the species status of taxa differentiated solely on 
chromosomal grounds have not been widely accepted 
(Sozen et al. 2006, Ivanitskaya et al. 2008). In accordance 
with the presentiy accepted view these taxa are herewith 
called "chromosomal forms". Alongside with taxonomic 
uncertainty the determination of conservation status of 
different mole rat taxa is further hampered by their exclu- 
sively subterranean lifestyle which makes it difficult to 
evaluate their population size. While the leucodon-super- 
species itself is categorised as Least Concern (Temple and 
Terry 2007), populations and habitats of many different 
European chromosomal forms are disappearing at an 
alarming rate, a phenomenon which has just recently been 
realized (Krystufek and Amori, 2008; Nemeth et al., 2009). 

In the course of researching mole rats of the Carpathian 
Basin, on both sides of the Hungarian-Serbian border a 
small and fragmented population of mole rats was identi- 
fied in 2008. Cytogenetic investigations have proven that 
this Kelebia-Suboticka pescara population belongs to one 
of the four endemic chromosomal forms of the Carpathian 
Basin (Nemeth et al. 2009), Nannospalax (leucodon) mon- 
tanosyrmiensis, described by Savic and Soldatovic (1974) 
as Spalax montanosyrmiensis. This form was previously 
known only from two neighbouring localities at the 
foothills of Fruska Gora in Serbia (called as Strazilovo- 
Cortanovci population). In this paper we provide data for 
population size, distribution area and threatening factors of 
this poorly known chromosomal form. 

Results and discussion 

Genetic identification 

Due to the methodological difficulties to culture very 
limited amount of blood cells for cytology we were able 
to karyotype a single individual only. The karyotype of 



a mole rat found in the vicinity of Kelebia was 2n = 54 
NF = 86 consisting of 2 pairs of metacentric autosomes, 
8 pairs of submetacentric autosomes, 5 pairs of subtelo- 
centric autosomes and 11 pairs of acrocentric auto- 
somes; the X chromosome is large and metacentric, the 
Y chromosome is medium sized and acrocentric. This 
karyotype is identical with that of previously reported 
from the vicinities of Strazilovo and Cortanovci in 
Serbia and described as montanosyrmiensis (Savic and 
Soldatovic 1974, Soldatovic and Savic 1983). 

As we were looking for evidence of the genetic differ- 
entiation of the montanosyrmiensis form and not the 
phylogeny of the superspecies only included endemic 
forms from the Carpathian Basin and the nominotypical 
form of the superspecies (Table 1); we also refrain from 
taxonomic conclusions. Based on the analysed 866 bp- 
long part of the cytochrome b gene sequences of three 
individuals of Vojvodina blind mole rat the degree of 
between-population sequence divergence with Kimura 
2-parameter distance was between 0.5 and 0.7%. This 
value is well below the average sister-species difference 
(which is more than 2.7%) in Rodentia (Bradley and 
Baker 2001). Thus the sequence analyses confirm the 
cytogenetic identity of the two populations. The sam- 
ples of the two populations of the Vojvodina blind mole 
rat form a discrete clade (Figure 1). This clade, with the 
approximate 10% p-distance, is well separated from other 
blind mole rat taxa inhabiting the Carpathian Basin and 
is also clearly distinct from specimens sampled at the 
type locality of the lesser blind mole rat Nannospalax 
[leucodon) leucodon. Sequences generated in the present 
study are deposited in GenBank under Accession Nos. 
JN656385-JN656390; tissue samples are stored in the 
Mammal Collection of the Hungarian Natural History 
Museum (HNHM) (Table 1). 

Recent distribution 

The Kelebia-Suboticka pescara population (Figure 2a) 
lives on both sides of the boundary between Serbia and 
Hungary. The Hungarian part (enclosed by the munici- 
pal boundaries of Kelebia) consists of 6 sub-populations 
isolated by roads, cultivated areas and wooded strips. 
These habitat fragments are small and deteriorated in 
quality (more than 50% of fallow land) covering a total 



Table 1 Blind mole rat specimens used in the genetic analyses 


Chromosomal form 


Settlement 


Country 


Voucher accession no. 


GenBank no. 


Montanosyrmiensis 


Kelebia 


Hungary 


HNHM 22789 


JN656386 


Montanosyrmiensis 


Strazilovo 


Serbia 


HNHM 23396 


JN656391 


Montanosyrmiensis 


Cortanovci 


Serbia 


HNHM 23397 


JN656392 


Leucodon 


Odessa 


Ukraine 


HNHM 23208 


JN656388 


Hungaricus 


Battonya 


Hungary 


HNHM 23001 


JN6S6387 


Transsylvanicus 


Hajdubagos 


Hungary 


HNHM 21838 


JN6S6385 
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Figure 1 Similarity dendrogram of Carpathian Basin blind mole rats and the nominotypical race from the type locality based on cyt b 
(866 bp) sequences. The dendrogram was constructed using UPGMA method based on similarity matrix calculated with "No. of differences" (the 
number of different base-pairs between two compared sequence). 



of 16 hectares. The estimated population size in 2010 
was approximately 90 individuals all together in the 6 
localities. The population on the Serbian side lives 
within the Suboticka pescara protected area. Although 
this contains 400 ha of temperate grassland (sand 
steppe), blind mole rats only occur on a 20-hectare 
stretch close to the international border. Despite the 
fact that this habitat (regenerated steppic grassland) is 
of relatively good quality the population numbers only 
about 60 individuals (Table 2). The reason could be that 
almost the entire area has been ploughed in 2003 which 
might have affected the population seriously. 

The Strazilovo-Cortanovci population (Figure 2b) is 
situated on the slopes of Fruska gora, almost 150 km 
south from the previous location. Clusters of mounds 
were found on pastures and meadows near Strazilovo 
and the population extends into the village of Corta- 
novci where one individual was caught in a kitchen gar- 
den. Blind mole-rats occur in a landscape mosaic of 
small woodlands, vineyards, orchards, pastures and 
meadows. Many of these cultivated habitat types have 
turned into spontaneously regenerating fallow land with 
elements of loess steppic vegetation. The actual habitat 
of the population is about 50 ha and the estimated 
population size is less than 100 individuals (Table 2). No 
part of the habitats of this population is protected. 

Landscape history 

The study area (10 x 10 km on the border of the settle- 
ments of Kelebia, Asotthalom and Suboticka) encloses all 
presently known fragments of the Kelebia-Suboticka pes- 
cara population. Following the categorisation of the mili- 
tary surveys types of open and closed dry grasslands, dry 
pastures and dry hayfields are regarded herewith as poten- 
tial mole rat habitats (see Topachevskii 1969, Savic and 
Nevo 1990, Horvath et al. 2007). 

Maps from the 1st military survey (Anon. 1785) shows 
a treeless, uninhabited region used for grazing livestock. 



In that era various open or more closed grassland com- 
munities covered 61% of the area. The 2nd military sur- 
vey (Anon. 1869) already shows profound changes. 
Spontaneously forming woodlands cover 15%, and sur- 
rounding the 67 homesteads cultivated land represents 
5% of surface cover; the potential blind mole-rat habi- 
tats decreased to 45%. A next survey in 1959 (Anon. 
1959) documented further increase in both woodland 
and cultivated area; drier grasslands formerly used as 
pastures were transformed by afforestation. Wooded 
areas totalled 35% in 1959 and their somewhat deceler- 
ating expansion is still ongoing into present days (Anon. 
2007). This affected much of the remaining grasslands 
the extent of which shrinked to 16% only and became 
heavily fragmented (Figure 3). 

Risk assessment 

In the Hungarian part of the Kelebia-Suboticka pescara 
population the major threat presently is the change in land 
management: tree plantations, small-scale farming and 
agro-industrial farming. Following the afforestation of the 
area in 2007, the extent of the habitat of one of the sub- 
populations was reduced by 50%. A small fragment (num- 
bering less than 20%) of the population falls into a Natura 
2000 site ("Deli Homokhdtsdg" SCI [Southern Sand Ridge 
SCI], code: HUKN20008), but the rest of the habitats are 
not protected by any means. In the Serbian part of the 
Kelebia-Suboticka pescara population the major threats 
are small-scale wood plantations. The whole area is pro- 
tected as part of the Suboticka pescara Protected Area, 
however, despite this status serious habitat destruction 
took place in 2003. 

In the Strazilovo-Cortanovci population the major threat 
is rapid habitat loss. Remaining grasslands and fallow lands 
are directly threatened. The field where the authors caught 
a specimen in 2009 was ploughed in 2010 wiping out all 
individuals on it. The subpopulations are isolated from 
each other and habitat fragments are under constant threat 
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(See figure on previous page.) 

Figure 2 Present distribution (shaded areas) of Vojvodina blind mole rat in the regions of a) Kelebia-Suboticka pescara and b) 
Strazilovo-Cortanovci. 



of transformation into cultivated land. Another obs- 
tacle to effective protection is that habitats are owned 
by a great number of proprietors, rendering long- 
term protection of the total population almost impos- 
sible. Presently, none of these areas are under any 
kind of protection. 

The montanosyrmiensis form is therefore proposed to 
be ranked as Critically Endangered according to the IUCN 
(2001) criteria. This assessment is based on the extent of 
occurrence which is estimated to be less than 100 km 2 ; 
area of occupancy is estimated to be less than 10 km . The 
populations are severely fragmented in no more than two 
locations and continuing decline is observed in area, 
extent and quality of habitats. 

Conclusions 

This paper presents a case study where a mammal 
within Europe can drift to the brink of extinction almost 
unnoticed as a result of the lack of information, unclear 
taxonomic status and unrecognised tasks in conserva- 
tion biology. 

Karyotype of the blind mole-rat populations flanking 
the boundary of Serbia and Hungary were found identi- 
cal to that of the montanosyrmiensis chromosomal 
form previously known only from Serbia (Soldatovic 



and Savic 1983). The mtDNA analysis confirmed 
that the two populations form a discrete clade which 
is highly separated from other blind mole-rat taxa 
found in the Carpathian Basin. According to the 
results of the most comprehensive molecular bio- 
logical research on blind mole rats based on 6 mito- 
chondrial genes, the montanosyrmiensis form is a 
well separated lineage that diverged from the closest 
taxon examined about 1.8 million years ago (Hadid 
et al. 2012). Therefore, montanosyrmiensis complies 
with the criteria of the Genetic Species Concept (Baker 
and Bradley 2006), as well as the ESU {Evolutionary 
Significant Unit) (Moritz 1994). 

Based on the surveys of all potential habitat patches 
in the presumed distribution area, the taxon montano- 
syrmiensis has only two extant populations that are 
150 km apart from each other and the total estimated 
number of individuals is less than 300. Even these two 
remaining populations are heavily fragmented and ex- 
tensive unsuitable vegetation types (forests, arable 
lands) between suitable habitat patches - often smaller 
than 3 hectares - possibly obstruct recent gene flow be- 
tween subpopulations. The status and future of these 
population fragments are precarious as 80% of the indivi- 
duals inhabits areas with no protection of any form, while 



Table 2 Recent populations of the Vojvodina blind mole rat 





Population 


Sub-population 


Extent (ha) 


Area of 
occupancy (ha) 


No. of individuals 
2009 2010 


Kelebia- 


Suboticka population 


Bacsborista 
Tanyahelyek 


5,5 
6,9 


2,1 
3,1 


12 
48 


19 
15 






Smuk-er (Tokleveles) 


9,3 


6,3 


55 


29 






Kelebiai halastavak mellett 


1,1 


0,5 


9 


10 






K6r6s-er 


2,2 


4,1 


24 


10 






Hatarszeglet 


2 


1,5 


12 


/ 






Suboticka pescara 


400 


20 


50 


60 






summa 


427 


37,6 


210 


150 


Strazilovo 


- Cortanovci population 


Strazilovo - fallow 


25,2 


18,3 


20-23 


10 






Strazilovo -meadow 


13,3 


10,2 


15-17 


16-18 






Strazilovo weekend houses 


4,2 


2,1 


12-14 


8-11 






vicinity of Cortanovci 


5,/ 


2,9 


10-12 


5-7 






close to the Fruska gora- Cortanovci road 


2,9 


1,6 


8 


6-8 






Cortanovci 


1,2 


1 


5-6 


5-6 






summa 


52,5 


36,1 


70-80 


50-60 






total population size 






280-290 


200-210 
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dry grassland | | quicksand 



Figure 3 Landscape history of the distribution area of the Kelebia-Suboticka pescara population based on four map sources 
(see details in text). Note the disappearence of prime blind mole rat habitats (dry grasslands and pastures) and the steady increase of 
unsuitable habitat types e.g. forests and cultivated fields. 



constantly being threatened by factors such as ploughing, 
afforestation and site development. As the landscape his- 
tory study revealed the Kelebia-Suboticka pescara popula- 
tion has lost more than 70% of its potential habitat in the 
last 200 years. 



Based on our present knowledge, the Vojvodina blind 
mole rat is one of the most seriously threatened, rarest 
mammal in Europe, the remaining population of which 
can be wiped out within years unless immediate conserva- 
tion action is taken. 
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Figure 4 Vojvodina blind mole rat (Nannospalax (leucodon) 
montanosyrmiensis) from Cortanovci, Serbia. 



Methods 

Distribution and population size 

Between 2008 and 2010 we systematically evaluated all 
known Hungarian and Serbian localities of Vojvodina 
blind mole rats (Figure 4) and also covered the remaining 
grasslands and fallow lands within its supposed distribu- 
tion area (Figure 5, Additional file 1). Only locations that 
were confirmed by the authors are listed below as "Recent 
distribution". Population size estimations were based on 
consistently monitoring the location of fresh mounds 
(Zuri and Terkel 1997). Since in all localities the clusters 
of mounds (one cluster is the visible sign of the under- 
ground activity of a single animal) were clearly separated 
from each other in this way it was possible to distinguish 
between mounds that belonged to different individuals. 
Counts were carried out in spring and in autumn months 
(during the peak activity periods of blind mole rats) in 
every locality at least in two consecutive years (Additional 
file 1). A population was defined as a set of sub-populations 
(inhabiting a more or less continuous, and homogeneous 
habitat) which were not separated from each other by 
impassable barriers (marshes or rivers). Gene flow between 
the two populations (i.e. the Kelebia-Suboticka pescara 
population and the Strazilovo-Cortanovci population) is 
quite improbable due to their distance and separating un- 
suitable habitats. 

Coordinates of all clusters of mounds were recorded 
by a hand-held GPS device. GPS records outlined the 
total extent of the distribution area. Spatial information, 
such as size and distribution of habitat patches, annual 
changes in the geographical position of mounds was 
tracked by means of Google Earth 5.1 and ArcGIS 9.0. 



Sampling 

As the Vojvodina blind mole rat is strictly protected 
both in Hungary and Serbia only 1-3 specimens were 
caught alive per sub-populations to prove the presence 
of the animals as allowed by research permits issued 
by the National Inspectorate for Environment, Nature 
Conservation and Water no. 14/05173-3/2006 and 14/ 
1840-3/2008. Altogether 16 animals were captured and 
handled in the field in accordance with guidelines 
approved by the American Society of Mammalogists 
(Gannon et al. 2007). The animals were caught by open- 
ing the tunnel system and capturing the animal trying to 
mend the damage (Nemeth et al. 2007). After samples 
were taken for tissue cultivation and DNA-analysis, indi- 
viduals were released at the site of capture straight into 
their own tunnels. For cytogenetic investigations, instead 
of using the lethal direct bone marrow preparation 
(widely used in rodents) we took blood samples aided by 
veterinarians, applying topical and systemic aenesthesia 
and 70% alcohol disinfection following the protocoll 
described by Sos et al. (2009). Lymphocyte cultures were 
used for chromosomal investigations. Metaphase arrest, 
hypotonic treatment and chromosome preparation were 
carried out according to standard cytogenetic techniques 
(Moorhead et al. 1960) with using RPMI 1640 medium 
containing 20% FCS, antimycotic, antibiotic solutions 
and the mixture of phytohaemagglutinin and pokeweed 
mitogen (1-1%). 

Sequence analysis 

The main objective of the molecular biological investiga- 
tion was to set the Kelebia-Suboticka pescara population 
and the Strazilovo-Cortanovci population in genetic 
comparison. Genomic DNA was isolated from fresh tis- 
sue samples with Qiagen DNeasy Blood and Tissue kit 
(Qiagen, California, USA) according to the extraction 
protocol. The 866 bp long sequences of cytochrome 
b was amplified with PCR using forward F-muarso 
(5'-ATGACATGAAAAATCATYGTTGT-3') and reverse 
R-muarso (5'-GAAATATCATTCKGGTTTAATRTG-3') 
primer pairs. The polymerase chain reaction was per- 
formed in total volume of 25 uL containing 30 ng tem- 
plate DNA, 1 uM of each oligonucleotide primer, 1.5 mM 
MgCl 2 , 160 uM dNTPs and 0.5U AmpliTaq DNA poly- 
merase (Applied Biosystems). PCR amplification was 
conducted in a DNA Engine Dyad (MJ Research) employ- 
ing 94°C for 2 min followed by 45 cycles of 94°C for 
45 sec, 48°C for 15 sec (ramp speed down to 48°C: lC/sec), 
60°C for 1 sec (ramp speed: 0.5C/sec up to 60°C) and 72°C 
for 2 min (ramp speed: lC/sec up to 72°C), and the final 
extension step was 72°C for 7 min. The PCR-product 
was checked on 1.6% agarose gel stained with ethidium 
bromide and cleaned with Qiagen Clean Up kit (Qiagen, 
California, USA). 
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Figure 5 Maps of investigated areas around the known populations (a and b) of Vojvodina blind mole rat. On the detailed maps yellow 
lines show international borders; red lines indicate the accessed routes; whereas yellow squares mark the throughoutly investigated potential 
habitat patches. 
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The A-tailing reaction was contained 7 uL of PCR prod- 
uct, 200 uM dATP, 1.5 mM MgC12, lx Taq buffer and 5U 
AmpliTaq DNA polymerase. The tailing profile was 94°C 
for 3 min and 70°C for 30 min. The A-tailed PCR product 
was ligated into pGEM Easy vector (Promega) and trans- 
formed into JM109 cells. Positive colonies were applied as a 
template for reamplification by PCR utilizing the T7 and 
SP6 primers. The products were checked on 1.6% agarose 
and then cleaned with SureClean (Bioline). The sequencing 
reaction was accomplished in GeneAmp PCR System 9700 
machine according to the thermal profile of 94°C for 4 min 
followed by 25 cycles of 94°C for 30 sec, 50°C for 15 sec 
and 60°C for 4 min. The 10 uL sequencing reaction con- 
sisted of 80 ng template, 2 uL of BigDye v3.1 Terminator, 
1 x buffer and 1 uM primer and then products were 
cleaned using the BigDye XTerminator Purification Kit 
Sequences for both directions were obtained using 
ABI3130 Genetic Analyser and were aligned and tested in 
MEGA v5 (Tamura et al. 2011). Prior to the analysis the 
aligned sequences were checked for the presence of stop 
codons by using ORF finder at NCBI (http://www.ncbi.nlm. 
nih.gov/gorf/gorf.html). Phylogenetic relation of the sam- 
ples was reconstructed using UPGMA method based on 
similarity matrix calculated with "no. of differences" (the 
number of different base-pairs between two compared se- 
quence) and "kimura 2-parameter" technique. P-distance 
was obtained by dividing the number of nucleotide differ- 
ences by the total number of nucleotides compared. 

Landscape history 

The landscape history could only be investigated on the 
two remaining populations of the Vojvodina blind mole 
rat but historical maps covering the last two centuries 
are available for the Kelebia-Suboticka pescara popula- 
tion only. These historical sources, which contain infor- 
mation on the changes of the landscape and of land use 
both in space and time, were the relevant digitized 
chart pages of military maps from the 1st (Anon. 1785) 
and 2nd (Anon. 1869) Military Mapping Survey of 
Austria-Hungary, and from the map series published by 
Honved Mapping Institute (Anon. 1959). Changes in 
land use during the past 50 years were assessed based 
on aerial photography (Anon. 2007). To track and in- 
terpret changes in landscape history ArcGIS 9.0 soft- 
ware was used. 



IUCN classifications and categories 

The Red List categories were assessed according to the 
2001 criteria (IUCN 2001). Habitats and threats were classi- 
fied (and terms used) according to the IUCN Habitats Clas- 
sification Scheme 3.0 and Threats Classification Scheme 
2.1, respectively (www.iucnredlist.org). 



Additional file 



Additional file 1: Surveyed potential sites within the presumed 
distribution area of the Vojvodina blind mole rat. 
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